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The voltage probe model is a model of incoherent scattering in quantum transport. Here we use 
this model to study the effect of spin-flip scattering on electrical conduction through a quantum 
dot with chaotic dynamics. The spin decay rate 7 is quantified by the correlation of spin-up and 
spin-down current fluctuations (spin- flip noise) . The resulting decoherence reduces the ability of the 
quantum dot to produce spin-entangled electron-hole pairs. For 7 greater than a critical value 7 C , the 
entanglement production rate vanishes identically. The statistical distribution P(7c) of the critical 
decay rate in an ensemble of chaotic quantum dots is calculated using the methods of random-matrix 
theory. For small y c this distribution is oc 7<T 1+/3,/2 , depending on the presence (fi = f) or absence 
(8 = 2) of time-reversal symmetry. To make contact with experimental observables, we derive a 
one-to-one relationship between the entanglement production rate and the spin-resolved shot noise, 
under the assumption that the density matrix is isotropic in the spin degrees of freedom. Unlike 
the Bell inequality, this relationship holds for both pure and mixed states. In the tunneling regime, 
the electron-hole pairs are entangled if and only if the correlator of parallel spin currents is at least 
twice larger than the correlator of antiparallel spin currents. 

PACS numbers: 73.23.-b, 03.67.Mn, 05.45.Mt, 73.63.Kv 
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I. INTRODUCTION 

The voltage probe model, introduced by Biittiker in 
the early days of mesoscopic physics, 1 gives a phenomeno- 
logical description of the loss of phase coherence in quan- 
tum transport. Electrons that enter the voltage probe 
are reinjected into the conductor with a random phase, 
so they can no longer contribute to quantum interference 
effects. Such a device is no substitute for a microscopic 
treatment of specific mechanisms of decoherence, but it 
serves a purpose in identifying model independent "uni- 
versal" features of the transition from coherent to inco- 
herent electrical conduction. 

In this work we introduce and analyze a novel appli- 
cation of the voltage probe model, to spin-resolved con- 
duction through a quantum dot. The voltage probe then 
serves a dual role: It randomizes the phase, as in the orig- 
inal spin-independent model, 1 but it also randomizes the 
spin. Two spin transport effects are examined: spin-flip 
noise and spin entanglement. The two effects are funda- 
mentally connected, in the sense that the degree to which 
spin-up and spin-down current fluctuations are correlated 
provides a measure of the degree of spin entanglement of 
electron-hole pairs exiting the quantum dot^ 

The geometry is sketched in Fig.^ The coupling of the 
electron spin to other degrees of freedom (nuclear spins, 
magnetic impurities, other electrons . . . ) is replaced by 
an artificial reservoir connected to the quantum dot via 
a tunnel barrier. The reservoir draws neither particles 
nor energy from the quantum dot£ Both the time aver- 
aged current and the time dependent fluctuations vanish, 
enforced by a fluctuating distribution function of the ar- 
tificial reservoir This phenomenological description of 
decoherence has found many applications in the context 
of (spin-independent) shot noise. (Recent references in- 



clude kUTI ISIOI . ) For alternative models of decoherence in 
that context, see Rcfs. 

In the context of spin-resolved conduction, the volt- 
age probe introduces two altogether different decay pro- 
cesses: spin flip and decoherence. These are character- 
ized in general by two independent decay times (denoted 
T\ and T2 < 2T±, respectively). In order to obtain two 
different time scales we could introduce, in addition to 
the spin-isotropic voltage probe, a pair of ferromagnetic 
voltage probes that randomize the phase without flipping 
the spin (pure dephasing) . Here we will restrict ourselves 
to the simplest model of a single voltage probe, corre- 
sponding to the limit T2 = 1T\ . This choice is motivated 
by the desire to have as few free parameters as possible 
in this exploratory study. The more general model will 
be needed to make contact with the existing microscopic 
theory for the spin decay times jiiji^iSiii 

The applications of the voltage probe model that we 
consider center around the concept of electron-hole en- 
tanglement. A voltage V applied over a single-channel 
conductor produces spin-entangled electron-hole pairs ^ 
The entanglement production rate is maximally eV/2h 
bits/second, for phase-coherent spin- independent scat- 
tering. Thermal fluctuations in the electron reservoirs^ 
as well as dephasing voltage fluctuations in the electro- 
magnetic environment^SiSi reduce the degree of entan- 
glement of the electron-hole pairs. Unlike other quantum 
interference effects, which decay smoothly to zero, the en- 
tanglement production rate vanishes identically beyond 
a critical temperature or beyond a critical decoherence 
rate. 

One goal of this investigation is to determine the prob- 
ability distribution of the critical decoherence rate in an 
ensemble of quantum dots with chaotic scattering. The 
fluctuations in the artificial reservoir reduce the entangle- 
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voltage probe 




FIG. 1: Illustration of the voltage probe model. A chaotic 
quantum dot is connected to a voltage source by two single- 
channel leads. Decoherence is introduced by means of a ficti- 
tious voltage probe, which conserves particle number within 
each energy range SE <C eV, on time scales St = h/SE. (The 
dashed line in the figure indicates the tunnel barrier that sep- 
arates the voltage probe from the quantum dot.) The random 
spin flips introduced by the voltage probe give a nonzero cor- 
relator of spin-up and spin-down currents (= spin-flip noise). 
The voltage probe also reduces the entanglement production 
by transforming the pure spin-singlet state of electron-hole 
pairs into a mixed Werner state. 



merit production by transforming the pure state of the 
electron-hole pair into a mixed state. For decoherence 
rates 7 > 7 C the density matrix of the electron-hole pairs 
becomes separable. The value of 7 C is sample specific, 
with a probability distribution P(^ c ) that we calculate 
using the methods of random-matrix theory^ 

The entanglement production is related to physical ob- 
servables via the spin-resolved shot noise. The correlator 
of spin-up and spin-down currents (spin-flip noise) is of 
particular interest, since it provides a direct measure of 
the spin relaxation time^ By assuming that the elastic 
scattering in the quantum dot is spin-independent (no 
spin-orbit interaction), we derive a one-to-one relation 
between the degree of entanglement (concurrence) of the 
electron-hole pairs and the spin-resolved shot noise. In 
the more general spin-dependent case such a relation ex- 
ists for pure states^iS^ through a Bell inequality, but not 
for mixed states^ The expressions for the concurrence C 
take a particularly simple form in the tunneling regime, 
where we find that C is nonzero if and only if the corre- 
lator of spin current fluctuations is at least twice larger 
for parallel spins than it is for antiparallel spins. 

We derive closed-form expressions for the ensemble av- 
eraged correlators in the regime of weak decoherence, 
both in the presence ((3 = 1) and absence ((3 = 2) of 
time-reversal symmetry. While the average spin-resolved 
current correlators are analytic in the decoherence rate 7 



around 7 = 0, the average concurrence (C) has a singular- 
ity at that point: a square-root singularity 1 — (C) oc ^7" 
for fj = 1, and a logarithmic singularity 1 — (C) oc 7I hi7| 
for (3 = 2. The singular effect of a small decoherence 
rate on the entanglement production also shows up in 
the probability distribution P(^ c ) of the critical decoher- 
ence rate: It does not vanish for 7 C — + 0, but instead has 
a large weight oc "f c 1+/3 / 2 . 

The outline of this paper is as follows. We start in Sec. 
ITT1 with a description of the system (quantum dot with 
voltage probe) and a formulation of the twofold question 
that we would like to answer (what is the entanglement 
production and how is it related to spin noise). A solution 
in general terms is presented in Sec. lIIII We begin in Sec. 
IIII Al bv simplifying the problem through the assumption 
of spin- independent scattering in the quantum dot. The 
concurrence of the electron-hole pairs is then given as 
a rational function of spin-resolved current correlators 
(Sec. IIII B|) . These correlators are expressed in terms of 
the scattering matrix elements of the quantum dot with 
voltage probe (Sec. IIII To evaluate these expressions 
an alternative formulation, in terms of a quantum dot 
without voltage probe but with an imaginary potential, is 
more convenient (Sec. lIIID|) . The connection to random- 
matrix theory is made in Sec. II VI By averaging over the 
random scattering matrices we obtain the nonanalytic 
7-dependencies mentioned above (Sees. Ivl and IVIfl . We 
conclude in Sec. IV11I 

II. FORMULATION OF THE PROBLEM 

We consider a quantum dot coupled to source and 
drain by single-channel point contacts. The voltage 
probe has channels, and is connected to the quantum 
dot by a barrier with a channel-independent tunnel prob- 
ability r^,. By taking the limit — > 0, — > 00 at fixed 
(dimensionless) conductance 7 = N^T^, we model spa- 
tially homogeneous decoherence with coherence time-Si 

Tcohoroncc = lim lim . (2.1) 

r^oJV^» 7A 

(We denote by A the mean spacing of spin-degenerate 
levels.) Since the mean dwell time in the quantum dot 
(without voltage probe) is Tdwcii = h/2A, one has 

7 ^T'dwell/' coherence* (^•^) 

The scattering matrix S of the whole system has di- 
mension (N^ + 2) x (N^ + 2). By convention the index 
n = 1 labels the source, the index n = 2 labels the drain, 
and the indices 3 < n < + 2 label the channels in the 
voltage probe. All of this refers to a single spin degree of 
freedom. Each channel is spin degenerate. As mentioned 
above, we assume that the scattering is spin independent. 
In particular, both the Zeeman energy and the spin-orbit 
coupling energy should be sufficiently small that spin ro- 
tation symmetry is not broken. The applied voltage V 
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between source and drain is assumed to be large com- 
pared to the temperature, but sufficiently small that the 
energy dependence of the scattering can be neglected. 

The energy range eV above the Fermi level is divided 
into small intervals SE -C eV. The voltage probe con- 
serves particle number and energy within each energy 
interval, on time scales St = h/SE. We write this re- 
quirement as I^E, t) = 0, where I^iE, t) is the electrical 
current through the voltage probe in the energy interval 
(E, E + SE), averaged over the time interval (t, t + 5t). 

Because the voltage probe does not couple different 
energy intervals, we may consider the entanglement pro- 
duction in each interval separately and sum over the in- 
tervals at the end of the calculation. In what follows we 
will refer to a single energy interval (without writing the 
energy argument explicitly). 

The density matrix p of the outgoing state in each en- 
ergy interval, traced over the degrees of freedom of the 
voltage probe, contains combinations of 0, 1, or 2 exci- 
tations in the spin degenerate channel of the source lead 
and the drain lead. Only the projection p e h onto a singly 
excited channel in the source as well as in the drain con- 
tributes to the entanglement production^ We denote by 
w = TiVp the weight of the projection, with V the op- 
erator that projects onto singly excited channels (so that 
wp c h — VpT^). The label "eh" stands for "electron-hole 
pair" , where "electron" refers to the single excited chan- 
nel in the drain and "hole" refers to the single nonexcited 
channel in the source. 

In the absence of decoherence /9 e h is pure (/9g h = p e h)- 
The voltage probe transforms p e h into a mixed state. Our 
aim is to calculate the loss of entanglement of p e h as a 
function of 7 and to relate it to the spin-resolved current 
correlators. 

III. GENERAL SOLUTION 

A. Simplification for spin-isotropic states 

The assumption of spin-independent scattering implies 
that the 4x4 density matrix p^ is invariant under the 
transformation (U <£> U)p h(U^ <8> U^) — p c h, for any 2x2 
unitary matrix U. As a consequence, p e h must be of the 
Werner formfSl 

Peh = i(l-Ol + C|*Beu)(*Bell|, -§<£<1, (3.1) 

with 1 the unit matrix and 

|*Bcn)=2- 1 / 2 (|U)-|IT)) (3.2) 

the Bell stated (The spin-up and spin-down arrows |, j. 
label the two eigenstates of the Pauli matrix a z .) 

The concurrence^ (degree of entanglement) of the 
Werner state is given by 

C = ^max(o,£-H. (3.3) 



The Werner state is separable for £ < 1/3. The entan- 
glement production rate (bits per second) in the energy 
range 8E under consideration is given byi& 

£ = ^^(i + iyi - ^), (3.4) 
T{x) = -xlog 2 x- (l-a:)log 2 (l-x). (3.5) 

The parameter £ that defines the Werner state 1)3. 1|) 
can be obtained from the spin-spin correlator 

TrpehCT z (g) a z = (3.6) 

In order to make contact with the voltage probe model 
we now relate this correlator to a spin-resolved current 
correlator, along the lines of Refs. l2 B24ll25l 



B. Solution in terms of current correlators 

We define N^^ft) as the number of electrons going 
out of the quantum dot in a time interval (t, t + St) 
through the source lead (X = S) or through the drain 
lead (X — D) with spin up (a =|) or with spin 
down (a =1). In terms of the current Ix a{t) one has 
N n? a (t) = -lD, a {t)St/e, Nffit) = 1 - IsAt)5t/e, with 
the convention that the current is positive if electrons 
enter the quantum dot. The spin-spin correlator (|3.6() is 
expressed by 

= ([Ng$(t) - N^(t)][N^(t) - N^\(t)}) 
([N^(t) + N^(t)][N^(t) + Ng%{t)\) 

= -(5t/e) 2 i<[J s>T (t)-J S4 (t)][7 D>T (t)-7 Ai (t)]), 

(3.7) 

w = (6t/e) 2 ([I s> i(t)+I s ,i(t)-2e/5t}[I Dt i(t)+I D , l (t)}), 

(3.8) 

where the brackets (• • • ) indicate an average over many 
measurements . 

The time dependent current Ix,a{t) — Ix.a + SIx, a (t) 
has time average Ix,a- The current fluctuations SIx, a (t) 
on the time scale St = h/SE have cross correlator 
(SIs,a{t)SlD,f3{t)) = (SE/h)P a p, with spectral density^ 



P a p= / dt(8Is, a (0)SI D ,e(t)). (3.9) 

J — oc 

The total spectral density of charge noise is given by 

/oo 
dt(5I s , a (0)8IsAt)) =-J2 P <*P- 
a.p -°° a,0 

(3.10) 

(The minus sign appears because J2p^s,/3 — 
— J2a^^D,/3, as a consequence of current conservation.) 
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Substitution into Eqs. I|3.7[l and (|3.8|l gives 

' {h/8E)(I s ,i - I s ,i)(Id,i - Id,i) 



e 2 SE w 



(3.11) 



(h/6E)(I s ,i + Is,i - 2eSE/h)(I DA + I D<i ) 

(3.12) 



e 2 5E 



Because of the spin isotropy, Ix,-\ = Ix.i, PfT = p lli 
and Pf j, = Pj.f . We denote by I > the total time aver- 
aged current from source to drain in the energy interval 
SE. Spin isotropy implies Is, a = \l and Id, a — 
Eqs. (|3~TT|) and [j3~T2j) then simplify to 



P 



1 ! 



P 



rr 



eI-\{h/5E)P 



w 



2h 
^SE 



1 p ■ 

2 r charge 
, f2 



(3.13) 



[el-i(V^)/ 2 - ^charge]- (3.14) 



C. Solution in terms of scattering matrix elements 

So far the analysis is general, not tied to a particular 
model of decoherence. Now we turn to the voltage probe 
model to express the average current and the current 
correlators in terms of the scattering matrix elements and 
the state of the reservoirs. For a recent exposition of this 
model we refer to Ref.0. The general expressions take the 
following form in the case of spin-independent scattering 
considered here: 



I = 

p n = 
p TT = 



2eSE ( T^T^z 

1 1 ,9 ~r 



h 

e 2 SE 

~2h~ { (N^-R^) 2 



e 2 5E 



Ns-Ri 



?12 + p n- 



(3.15) 

(3.16) 
(3.17) 



and the correlators of intrinsic current fluctuations^ 

Qnm ^ (J^nn'&nm 1 ^iiti'^™') 

n' ,m' — 1 

X &mn' ^rarn' ^mn') fn' (1 /m')}(3.20) 

N^+2 

Qn<f> ^ Qnmi Q(jxf> — ^ Qnm- (3-21) 



m— 3 



n,m— 3 



The state incident on the quantum dot from the reser- 
voirs is fully characterized by the mean occupation num- 
ber /„, given by 



1 if n = 1, 

_ . if n = 2, 

In - \ T 



N<t> — Rd 



if 3 < n < + 2. 



(3.22) 



For the source and drain reservoirs this is a state of ther- 
mal equilibrium at zero temperature. For the fictitious 
reservoir this is the nonequilibrium state 

P,p= II (/n4|0)(0K + (l-/ n )|0)(0|), (3.23) 

n— 3 

with aj, the operator that excites the n-th mode in the 
voltage probe. These are all Gaussian states, meaning 
that averages of powers of a n and can be constructed 
out of the second moment (a^an) = f n by the rule of 
Gaussian averages. 



D. Reformulation in terms of imaginary potential 
model 



The model of a quantum dot with voltage probe can 
be reformulated in terms of a quantum dot without volt- 
age probe but with an imaginary potential^ This refor- 
mulation simplifies the expressions for the entanglement 
production, so we will carry it out here. 

The unitarity of S makes it possible to eliminate from 
the expressions in Sec. IIII 01 all matrix elements that in- 
volve the voltage probe. Only the four matrix elements 
S nm , n,m G {1,2}, involving the source and drain re- 
main. This sub-matrix of S forms the sub-unitary matrix 



We have defined the transmission and reflection 
probabilities 1 



(Sn 


Sl 2 \ 




3 


[s 21 


S22 J 


"I 





(3.24) 



Pn — >m — I^Vnnl , R^ — ^ ^ I'S'nml , 



n,m— 3 



N^+2 



(3.18) 



P(p— >m — ^ \Srmi\: Pji^di — ^ \Smn 1^,3.19) 



n— 3 



m— 3 



As derived in Ref. |2y, the matrix s corresponds to the 
scattering matrix of the quantum dot without the voltage 
probe, but with a spatially uniform imaginary potential 
— ijA/Air. The coefficients t,t',r,r' are the transmission 
and reflection amplitudes of the quantum dot with the 
imaginary potential. 
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After performing this elimination, the expressions 
(|3.20hl3~2TT)l for the correlators Q nm that we need take 
the form 

0i2 = -\Ci--U)s 11 s; 1 -us 12 s; 2 \ 2 , (3.25) 

Qll - [/0(l-|^12| 2 ) + (l-/0)|^ll| 2 ] 

x [(1-/ )(1-|5„| 2 ) + / |5 12 | 2 ] ; (3.26) 
= 2(Q 12 + Q U ), (3.27) 
Qi4> = Q24, — ~^Q<p<p, (3.28) 

with mean occupation number 

= T w = 1- (at a ) n 
N^-R^ 2 - Tr sst ' 

We also define the quantity 

7 _ 1 - {ss r )n 



(3.29) 



Ns-Rd 



2 - Tr sst 



(3.30) 



which equals in the presence of time-reversal symme- 
try (when T, wm = T m ^ n ) — but is different in general. 

The expressions (|3.15fl3~T7jl for the mean current / and 
the correlators P a p simplify to 

/= ^[Ua-\S22\ 2 ) + (l-U)\S 21 \ 2 )}, (3.31) 



p 



e 2 SE „ e 2 SE 



11 



?X2 = 



2h 



Q<t>d,[U{i - U) - 5]. 

(3.32) 



Some more algebra shows that 



el - \{h/8E)P = 2{e 2 8E/h)Q lll (3.33a) 
Q H = 2/^(1 -/ )(1-Det sst). (3.33b) 

Substitution of Eqs. and (|3~33|) into Eqs. <PTT3|) 

and 13.14|l finally gives compact expressions for the 
Werner parameter £ and the weight w of the electron- 
hole pair: 



£ = 



(3.34) 
(3.35) 



X + Y 
w = 2(X + Y), 

X = / (l-/ )[2/ (l-/ ) + l] 

x(l-Dctsst), (3.36) 

Y = \rt* - U( SS ^ 12 \ 2 . (3.37) 

The spin-resolved current correlators 1)3. 32[) are expressed 
similarly by 



Pu = P 
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h ' h K2 
Z = / (1-/ )(1-Detsst) 



V —^[\X-Z), (3.38) 
(3.39) 



Let us check that we recover the known result for the 
entanglement production in the absence of decoherence. 



In that case s is a unitary matrix sq, so X, Z — > and 
Y — > |r tol 2 — independent of (The label indi- 
cates zero decoherence rate.) Hence £ = 1 (maximally 
entangled electron-hole pairs) and 



w = 2g a (l -g ), 



(3.40) 



with go = \to\ 2 the phase coherent conductance of the 
quantum dot in units of 2e 2 /h. The total entanglement 
production rate (integrated over all energies) becomes 



So = (eV/h)w = (2eV/% (l -50), 



(3.41) 



in agreement with Ref. 0- Furthermore, we verify that 
in this case P^ — (no spin- flip scattering without 
the voltage probe), while P^ = ~(e 2 8E/h)g (l — go) 
is given by the shot noise formula for spin-independent 
scattering. 32 ! 33 



IV. RANDOM-MATRIX THEORY 
A. Distribution of scattering matrices 

The expressions of the previous section refer to a single 
quantum dot. We now consider an ensemble of quantum 
dots, generated by small variations in shape or Fermi en- 
ergy. For chaotic scattering the ensemble of scattering 
matrices is described by random-matrix theory, charac- 
terized by the symmetry index [3 = 1 in the presence of 
time-reversal symmetry and (3 = 2 if time-reversal sym- 
metry is broken by a magnetic field. 22 (The magnetic 
field should be sufficiently weak that the Zeeman energy 
does not lift the spin degeneracy.) Since we assume that 
spin-orbit coupling is not strong enough to break the spin 
rotation symmetry, the case (3 — 4 of symplectic symme- 
try does not appear^ 

In the absence of decoherence, s is unitary and its dis- 
tribution is the circular ensemble. With decoherence, s is 
sub-unitary. Its distribution P(s) was calculated in Ref. 
I2II It is given in terms of the polar decomposition 







y/T=l£ 



(4.1) 



with unitary matrices v! = u T if (3 = 1 and v! inde- 
pendent of u if f3 = 2. These matrices are uniformly 
distributed in the unitary group. The real numbers 
Tx,T2 € [0, 1] are the two eigenvalues of 1 — sst. Their 
distribution Pa (t\ , T2 ) is given as a function of 7 by Eq. 
(17) of Ref. 26. (It is a rather lengthy expression, so we 
do not repeat it here.) 

We parameterize the 2x2 unitary matrix u by 



sin a 



sina 



,-Wl-!Q2 



cos a 



(4.2) 



and similarly for u' . The angles a\, a%, 03 are uniformly 
distributed in the interval (0, 2tt), while the angle a € 
(0,7r/2) has distribution P(a) — sin 2a. 
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In this parameterization the occupation numbers with a' = a if j3 = 1. The quantities X, Y, Z that 



(|3~2"9"|) and i|3~3U|l arc 



determine £, w, P^, P^ become 



5 + | cos 2a 



, ti - r 2 

Tl + T 2 ' 



> , i cos 2a 



Tl + T2 



(4.3) 
(4.4) 



X = U(1-U)[2U(1-U) + 1}(t 1 +T2-t 1 t 2 ), 



Y 



Vl — ti sin a cos a' — e l * yT — t 2 cos a sin a') (e l *-\/l — Ti cos a cos a' + e ** %/l — t 2 sin a sin a') 



+ \U{t\ - T 2 )sin2a , 
Z = f<f,(l - /^)(ti +t 2 -t x t 2 ). 



(4.5) 

) 

(4.6) 
(4.7) 



The phase $ = a 2 +a^ is uniformly distributed in (0, 27r), 
regardless of the value of (3. 



B. Weak decoherence 



In the regime 7 < 1 of weak decoherence the expres- 
sions simplify considerably. The distribution P^(t\,ti) 
is then given by the Laguerre ensemble^S^ 



Tl -T 2 | 



(TlT 2 )2/3+2' 

(4.8) 

with ci = 1/48 and c 2 = 1/24. Since t 15 t 2 < 7 <C 1, we 
may expand 



X 



U(i ~ /*)[2/*(i - /*) + i](n + t 2 ) + o(Tf ), 



(4.9) 

^ = 90(1 - flb)(l - Tl - T 2 ) + (Tl - T 2 )(/ - i) 

x [(.90 - |) cos 2a + i cos 2a'] + 0(rf ), (4.10) 

Z = / (l-/ )(T 1+ T 2 )+0(Tf). (4.11) 



The phase coherent conductance go = |io| 2 is given in 
terms of the angles a, a', <& by 

.g = I — I cos 2a cos 2a' — | sin 2a sin 2a' cos 2$. (4.12) 

It is independent of n and r 2 , with distribution 36 



^(<?o) = 5^ 1+/?/2 : < .go < 1 



(4.13) 



V. ENSEMBLE AVERAGES 

Averages over the ensemble of chaotic cavities require 
a four- fold integration for (3 — 1 (when a' = a), 

1 r l ,2tt ,tt/2 

dridr 2 Pi(ri,T 2 ) / — / sin 2a da • ■ 

(5.1) 



27r jo 



and a five- fold integration for (3 = 2, 



•••/2 = 



,tt/2 W2 

x / / sin 2a sin 2a' dada' 
J Jo 



(5.2) 



Results are plotted in Figs. [3 and |3 In Fig. [21 we see 
that the correlator Pjf of parallel spin currents (lower 
curves) is reduced in absolute value by the voltage probe 
— in contrast to the spin-flip noise Pfj. (upper curves), 
which is increased in absolute value. For large 7 all cor- 
relators tend to the same limit, 



lim P aa > = — , 

7^00 16 h 



(5.3) 



regardless of the presence or absence of time-reversal 
symmetry. In Fig. [3] we see how the decoherence intro- 
duced by the voltage probe reduces both the entangle- 
ment per electron-hole pair (quantified by the concur- 
rence C), as well as the total entanglement production 
rate £. 

In the limit of weak decoherence, the averages can be 
calculated in closed form using the formulas from Sec. 
HVBl For the spin-resolved current correlators we find, 
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to order j 2 : 
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FIG. 2: Dependence on the dimensionless decoherence rate 
7 = 2rdwcii/Tcohcrcncc of the ensemble averaged spin-resolved 
current correlators P-\[ and P TT , both in the presence (/3 = 1) 
and absence (/3 = 2) of time-reversal symmetry. The solid 
and dashed curves are computed by averaging Eq. 13.3811 with 
the random-matrix distributions, according to Eqs. 115. Ill and 
15.211 . The dotted lines are the weak-decoherence asymptotes 
15.4H5.71 . For strong decoherence all curves tend to the value 
-JE^V/h. 




FIG. 3: Dependence on 7 of the average concurrence C and 
entanglement production rate £ . The solid and dashed curves 
are computed by averaging Eqs. l|3.3H3.4t and l|3. 34143.3"^ . 
The dotted lines are the weak-decoherence asymptotes 15.91 - 
The asymptote for (C) 2 converges poorly, because the 
next term of order 7 in Eq. 15,101 is not much smaller than 
the term of order 7 In 7. 
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(5.4) 
(5.5) 
(5.6) 

(5.7) 



(We have replaced 8E by eV, to obtain the total inte- 
grated contributions.) 

The average Werner parameter, 



(Op = 1 - 



x 



X + Y 



(5.8) 



is nonanalytic in 7 around 7 = 0, because (X/go)p di- 
verges. Since P(go) oc g a 1+/3 / 2 ; c f. Eq. I|4.13|l . the average 
has a square-root singularity for [3 = 1 and a logarithmic 
singularity for (3 — 2. The average concurrence has the 
same singularity, in view of Eq. (|3.3|) . To leading order 
in 7 we find 



(£>i 
(O2 



0.75 7 1/2 

13 1 
— 7 In — 

42 ' 7 



(C)i 

(C>2 



1.13 7 1/2 , (5.9) 
|7lni(5.10) 



The ensemble averaged weight (w)p = 2(X + Y)p of 
the electron- hole pairs is analytic in -) , 

, . 4 11 , . 1 62 , 

(w)l = 1 7, (w)n = 1 7. (5.11) 

w 15 30 v 1 3 189' 1 ' 

The average entanglement production is given, to leading 
order in 7, by 



eV 
~~h 



2(X + Y) I3 
4 11 



( ^ll5 + 30 7 - 
w \3 189 ' 631n2 




(5.12) 



(5.13) 



(We have again replaced SE by eV for the total entan- 
glement production.) 



VI. CRITICAL DECOHERENCE RATE 

For each quantum dot in the ensemble, the entangle- 
ment production rate £ vanishes identically for 7 greater 
than a certain value 7 C at which the Werner parameter £ 
has dropped to 1/3. For 7 slightly less than j c we may 
expand £ = 1/3 + 0("f c — 7). In view of Eqs. I|3.3|l and 
(|3.4|) the entanglement production rate has a logarithmic 
singularity at the critical point, 



E cx (7 C - 7) 2 In (7 C - 7) \ if 7 | 7c 



(6.1) 
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FIG. 4: Probability distribution of the critical decoherence 
rate 7 C , beyond which the entanglement production vanishes. 
The solid and dashed curves are computed from Eq. i10.2> . 
The dotted lines are the weak-decoherence asymptotes 16.41 : 
P(jc) oc 7c " 1+/3/2 for 7c < 1. 



This is a generic feature of the loss of entanglement by the 
transition to a mixed state, cf. the logarithmic singular- 
ity in the temperature dependence of the entanglement 
production found in Ref. 

The statistical distribution Pp(j c ) of the critical deco- 
herence rate in the ensemble of chaotic quantum dots is 
defined by 



Pfii-Yc) 



(6.2) 



with Q(x) the unit step function (Q(x) = 1 if x > and 
Q(x) — if x < 0). The result of a numerical evaluation 
of Eq. (|6.2|l is plotted in Fig.QJ The ensemble average is 



j 0.954 if 3 



1, 
2. 



(6.3) 



Since 7 = 2rd W eii/T C oherenco, the critical decoherence rate 
of a typical sample in the ensemble of chaotic quantum 
dots is of the order of the inverse of the mean dwell time. 
Although the mean of the distributions for 8 = 1 and 
3 = 2 is almost the same, their shape is entirely different, 
cf. Fig.H 

The full probability distribution shows that sample-to- 
sample fluctuations are large, with a substantial weight 
for 7 C -c 1. For small 7 C the distribution Ppdc) has the 



same limiting behavior oc 7, 



-1+/3/2 



as the conductance go 



[cf. Eq. Ij4.13|l ]. More precisely, as derived in Add. lAl 



,. j, f \ J 0.085 7t 
^0^ = 113/42 



-1/2 



if p = 1, 
if 3 = 2. 



(6.4) 



A. Strength and weakness of the voltage probe 
model 



We have shown how the voltage probe model of shot 
nois oMi 5 can be used to study spin relaxation and deco- 
herence in electrical conduction through a quantum dot. 
The strength of this approach to spin transport is that 
it is nonperturbative in the dimensionless conductance 
go, permitting a solution for go of order unity using the 
methods of random-matrix theory. It is therefore com- 
plementary to existing semiclassical approaches to spin 
noise, 3 '' which require go 3> 1. 

The weakness of the voltage probe model is that it 
is phenomenological, not directly related to any specific 
mechanism for decoherence. We have examined here the 
simplest implementation with a single voltage probe, cor- 
responding to a single decay rate 7. The dominant deco- 
herence mechanism of electron spins in a quantum dot, 
hyperfine coupling to nuclear spins^ has a much shorter 
(ensemble averaged) decoherence time T2 than the spin- 
flip time T\. Pure dephasing (decoherence without spin 
flips) can be included into the model by means of ferro- 
magnetic voltage probes. This is one extension that we 
leave for future investigation. 

Another extension is to include spin-orbit scattering 
(symmetry index 3 = 4). We surmise that the result 

P (7 C ) oc 7c 1+l3 ^ 2 for the distribution of the critical de- 
coherence rate in the weak decoherence regime, derived 
here for the case 8 = 1, 2 without spin-orbit scattering, 
holds for 8 — 4 as well — but this still needs to be 
demonstrated. 



B. Entanglement detection for spin-isotropic states 

By restricting ourselves to a system without a prefer- 
ential basis in spin space, we have derived in Sec. IIII Bl 
a one-to-one relation between the entanglement produc- 
tion and the spin-resolved shot noise. This relation goes 
beyond the voltage probe model, so we discuss it here in 
more general terms. 

The basic assumption is that the conduction electrons 
have no preferential quantization axis for the spin. This 
socalled SU(2) invariance means that the full density ma- 
trix p is invariant under the simultaneous rotation of each 
electron spin by any 2x2 unitary matrix U: 

U ®U ®U ■ ■ ■f&UpU't ® •••[/ t ®U ] =p. (7.1) 

The 4x4 matrix p c h, obtained from p by projecting onto 
a single excited channel in the source as well as in the 
drain, has the same invariance property: 

U <E) Up ch U f ® [7 f = peh. (7.2) 

As explained in Sec. IIII Al the concurrence of the 
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C = - max { 0, -Tr p ch a z 



(7.3) 



electron-hole pairs then follows directly from 

1 

z_ 3 

Here we have excluded a spontaneous breaking of the 
spin-rotation symmetry (no ferromagnetic order). The 
more general case has been considered in the context of 
the isotropic Heisenberg models 

The concurrence in this spin-isotropic case is related 
to the spin-resolved shot noise by Eqs. I|3.3|l and l|3.13|l . 
The entanglement production rate £ follows according to 
Eq. (|3.4|l from C and a weight factor w, given in terms 
of the shot noise by Eq. (|3.14|) . To detect the spin en- 
tanglement one thus needs to measure the correlator of 
parallel and anti-parallel spin currents. This is in essence 
a form of "quantum state tomography" , simplified by the 
fact that an SU(2) invariant mixed state of two qubits is 
described by a single parameter (the Werner parameter 
£). The isotropy assumption does away with the need to 
compare correlators in different bases, as required for the 
Bell inequality 24,25 or for quantum state tomography of 
an arbitrary density matrixes 

In closing, we mention the remarkable simplification of 
the general expressions of Sec. lIIIBl if the (dimensionless) 
conductance go <C 1 (tunneling regime). Then the shot 
noise is Poissonian, hence P c har g e = el = — 2(Pff + PfJ.)- 
Moreover, the term quadratic in / is smaller than the 
term linear in / by a factor go , so it may be neglected. 
Instead of Eqs. (|3.13|l and l|3.14|l we thus have 



P 
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p 



w 



z5E 



(7.4) 



P TT + P u 

The expressions ()3.3Jl and l|3.4(l for the concurrence and 
entanglement production simplify to 



C = — max < 0, 



Pi 



TT 



P 



n 



£ 



P TT H 
£=max{0,|P TT | 



It 

e 



P 



11 



(7.5) 
(7.6) 



(7.7) 



We thus arrive at the conclusion that the electron-hole 
pairs produced by a tunnel barrier in a single-channel 
conductor with spin-independent scattering are entan- 
gled if and only if \P^\ > 2|P|-jJ, that is to say, if and 
only if the correlator of parallel spin currents is at least 
twice as large as the correlator of antiparallel spin cur- 
rents. We hope that this simple entanglement criterion 
will motivate further experimental efforts in the detection 
of spin noise. 
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APPENDIX A: DERIVATION OF EQ. (ET^ 



expressions of Sec. lIVBl for the weak decoherence regime. 

If 7 <C 1 the criticality condition £ = 1/3 is equivalent 
to <7o(l — go) = Ql, with the definition 

Q = (ri+r 2 )[U(l- U)(m~ U) + l)+9o(l~9o)] 
- (n - t 2 ){U - |) [(So - \) cos 2a + \ cos 2a']. 

(Al) 

The Laguerre distribution l|4.8|l of the rescaled variables 
Ti = Ti /7 is independent of 7 in the limit 7 — > (when f% 
ranges from to 00). Substitution into Eq. (|6.2() gives 



n , , / ? , f 9oO--9o) 
Pp{lc) = { 8 ( 7, 



(A2) 







We first consider the case j3 = 1. Then a' = a, so go 
and Q simplify to 

50 = (sin 2a sin $) 2 , (A3) 
Q = (ri+T 2 )[/ 2 (l-/ ) 2 + i/ (l-/ ) + « ?o (l- 5o )] 
- (n - f 2 )(U - i)gocos2a. (A4) 

The average over $ contributes predominantly near $ = 
and <I> = 7r, with the result 



lim Pi (7 C ) = —\ — 
7c^o 2tt V 7c 



, (r 1 +f 2 ) 1 / 2 [/|(l-/ ) 2 + ^(l-/,)] 
sin 2a 



1/2' 



(A5) 



The remaining average over and a gives simply a nu- 
merical coefficient, resulting in Eq. I|6.4|l . 

Turning now to the case f3 = 2, we first observe that 
the limit j c — > contains equal contributions from go 
near and 1. Hence Eq. (|A2|) simplifies to 



(A6) 



lim P 2 ( 7c )=2(Q%o)) 2 . 

7c ^0 



To reach go — we need a = a' and $ = or n. Expand- 
ing around a = of and $ = 0, we have to second order 
go = < E )2 sin 2 2a + (a — a') 2 . There is a similar expansion 
around $ = ir. Using the identity S(a 2 + b 2 ) = tt 5 (a) 6(b) 
we thus arrive at 



S(9o) 



-6(a-a')[5($) + 6($-Tr)]. (A7) 



sin 2a 

Substitution into Eq. i|A6(l gives the limiting value 
nmP 2 (7c) = ((fi+r 2 )[2f 2 (l-U) 2 + U(l-U)]) 2 
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42' 

as stated in Eq. I|6.4|) . 



(A8) 



We wish to evaluate the distribution Pp^c) of the crit- 
ical decoherence rate in the limit 7 C — > 0. We can use the 
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